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^ ■ Classical novae are fascinating stellar events, at the crossroads of astrophysics, nuclear 

physics and cosmochemistry. In this review, we outline the history of nova modeling with 
^ . special emphasis on recent advances and perspectives in multidimensional simulations. 
OO I Among the topics that are covered, we analyze the interplay between nova outbursts 
and the Galactic chemical abundances, the synthesis of radioactive nuclei of interest for 
■ gamma-ray astronomy, such as ''Li, ^^Na or ^^Al, and the recent discovery of presolar 
Q . meteoritic grains, likely condensed in nova shells. 
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^ . 1. Nuclear ashes: Classical novae and Galactic nucleosynthesis 
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Classical novae are close binary systems consisting of a white dwarf, and a large main 
sequence (or a more evolved) star. The companion overfills its Roche lobe and matter 
^ I flows through the inner Lagrangian point, leading to the formation of an accretion disk 
around the compact star. A fraction of this (H-rich) matter ultimately ends up on top of 
the white dwarf, where it is gradually compressed up to the point when ignition conditions 
to drive a thermonuclear runaway (hereafter, TNR) are reached. 
5^ I The thermonuclear origin of nova outbursts was first theorized by Schatzman |3ti|37j . 
Modern multiwavelength observations and numerical simulations (pioneered by the early 
hydro models of Starrfield et al. [10]) have drawn a basic picture, usually referred to as 
the thermonuclear runaway model. Since then, several groups have attempted to improve 
our understanding of these dramatic stellar events, including state-of-the-art nova nucle- 
osynthesis studies with spherically symmetric (or 1-D) hydro codes (see |2(j|21|43j . and 
references therein) or preliminary multi-D approaches jl3|14|24l25j . 

Nuclear physics plays a crucial role in the course of the explosion. As material from the 
accretion disk piles up on top of the (CO or ONe) white dwarf, the first nuclear reactions 
take place. This follows a rise in temperature since degenerate conditions unable the star 
to readjust the hydrostatic equilibrium by an envelope expansion and, as a result, a TNR 
ensues. The triggering reaction is ^^C(p,7), which initiates the 'cold' CNO cycle. At very 
early stages of the explosion, the main nuclear activity is driven by ^^C(p, 7)^^N(/5"'")^^C(p, 
7)^^N. But as the temperature rises, the characteristic time for proton capture reactions 
on ^^N becomes shorther than its (3~^ decay time, initiating the 'hot' CNO cycle. This 
is accompanyied by proton capture reactions onto ^^N, leading to ^^O, as well as by 
^^0(p,7)^''F (near peak temperature, Tpeak)- At this stage, the envelope exhibits the 
presence of significant amounts of ^^N, ^^'^^O and ^^F. Indeed, it is the decay of these 
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short-lived, (3^ unstable nuclei, carried away from the hot envelope's base to the outer, 
cooler layers by convective transport, what powers the ejection phase ^U]: their sudden 
release of energy, few minutes after the peak of the explosion, increases the entropy and 
temperature of the material. As a result, the electron degeneracy is lifted and an overall 
expansion sets in, driving the ejection of most (if not all) of the accreted envelope into the 
interstellar medium (indeed, the TNR is halted by the expansion of the envelope rather 
than by fuel consumption). The ejected layers enclose the history of multiple nuclear 
processes that modified its chemical composition during the course of the TNR: they are 
characterized by huge amounts of the daughther nuclei ^^N, ^^O, and to some extent, 
^^C, which have been claimed to represent the fingerprints of classical nova outbursts in 
the overall Galactic history f |26|2H4nj . and references therein). Moreover, the envelope 
is also enriched in other species (depending on the nova type, the mass-accretion rate or 
the white dwarf initial luminosity), that may contribute to the Galactic abundances to a 
lesser extent. This includes nuclei such as ^Li, ^^F, or ^^Al. 

In order to match the energetics, peak luminosities, and associated nucleosynthesis of 
the so-called fast novae, the modeling of the explosion requires mixing between the solar- 
like material transferred from the companion and the outermost layers (CO- or ONe-rich) 
of the underlying white dwarf. In fact, at the typical temperatures expected during 
the course of the TNR, the amount of leakage from the CNO cycle is very limited |21j . 
and hence, the observed abundances of elements ranging from Ne to Ca (significantly 
overproduced with respect to solar proportions in some novae) cannot be explained in a 
natural way in terms of nuclear processes. Indeed, the quest for a self-consistent mixing 
mechanism has become the Holy Grail of nova modeling: several mechanisms have been 
proposed so far, including diffusion induced convection |3Ull9j . shear mixing [29 28j, con- 
vective overshoot induced flame propagation ^3], convection induced shear mixing |27j . 
or more recently, mixing by gravity wave breaking on white dwarf surfaces |33|lj . But the 
final word has not yet been said... 

Among the species synthesized during classical nova outbursts, several radioactive nuclei 
have deserved special attention, in particular those associated with the predicted gamma- 
ray output from novae. Whereas ^^N and ^^F may be responsible for a prompt 7-ray 
emission at and below 511 keV, ''Be ("^Li) and ^^Na |15|17j . which decay much later, 
may power line emission at 478 and 1275 keV, respectively. ^^Al is another interesting 
radioactive isotope that can be synthesized during nova outbursts, although due to its long 
lifetime only its cumulative emission can be observed. A detailed account of the predicted 
gamma-ray signatures of classical nova outbursts (including detectability distances) can 
be found in ^H]. Hereafter, we will concentrate on the most relevant aspects of nova 
modeling, pointing out where do we stand and what is missing in our understanding of 
the overall phenomenon. In particular, we will focus on the synthesis paths of ^Li, ^^Na 
and ^^Al, the associated nucleosynthesis and observational constraints. 

1.1. ^Li 

Recent hydrodynamic simulations |16|21j have confirmed the feasibility of the Be- 
transport mechanism |5| as the key for ^Li production in nova outbursts: the process 
is initiated by the synthesis of '^Be through ^He(a, 7) ''Be, which is ultimately transformed 
into ^Li (r ~ 77 days) by means of an electron capture, with the emission of a charac- 
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teristic 478 keV 7-ray photon. Huge ^Li overproduction factors (i.e., ~ 900) with respect 
to solar values have been obtained, in particular for novae hosting CO white dwarf cores. 
These studies stressed the critical role played by the quasi-equilibrium between ''Be and 
^B, driven by efficient photodisintegration reactions on ^B, on the survival of ^Be around 
peak temperatures, and confirmed novae as likely ^Li factories [IT] (earlier hydrodynamic 
simulations assumed, however, envelopes in-place, thus neglecting the possible impact of 
the initial stages of the TNR and the onset of convective transport on ''Li production). 
Moreover, they refuted the conclusions based on simple 1- and/or 2- zone models ^ that 
stressed the key role played by ^B(p,7) in breaking the quasi-equilibrium between ^Be 
and ^B (and hence, leading to ''Be destruction instead). No nuclear uncertainties in the 
domain of nova temperatures significantly affect ''Li synthesis |16j . 

It is worth noting that the potential contribution of classical novae to the Galactic ''Li 
content turns out to be small (i.e., less than 15%, according to |16|21j ). However, a nova 
contribution seems to be required to match the ''Li content in realistic calculations of 
Galactic chemical evolution |32j . 

1.2. 22Na 

The role of ^^Na for diagnosis of nova outbursts was first outlined in the seminal work 
of Clayton & Hoyle [6^ . This isotope decays into a short-lived excited state of ^^Ne, which 
de-excites to its ground state by emitting a 7-ray photon of 1.275 MeV. 

The synthesis of ^^Na in novae proceeds through two alternative reaction paths. In the 
Ne-enriched envelopes of ONe novae it takes place through ^°Ne(p,7)^^Na, followed 
either by another proton capture and then, a /5+-decay into ^^Na, ^^Na(p,7)^^Mg(/?"'")^^Na, 
or decaying first into ^^Ne before another proton capture ensues, ^"'^Na(/5"'")^^Ne(p,7)^^Na. 
The main destruction channel at nova temperatures is ^^Na(p,7)^^Mg. 

The nuclear uncertainties associated with the synthesis of ^^Na in novae [22] have 
been recently reduced due to the first direct measurement of the ^^Na(p,7) rate with 
the DRAGON recoil separator at TRIUMF |3|8j . and to indirect determinations of the 
^^Na(p,7) rate carried out with the Gammasphere at the Argonne National Lab ^20] . 

1.3. 26 Al 

2^A1 was discovered in the interstellar medium by the HEAO-3 satellite, through the 
detection of the 1809 keV 7-ray line. This characteristic 7-ray feature is produced by the 
decay (r = 1.04 Myr) of the ^^Al ground state to the first excited state of ^^Mg, which 
in turn de-excites to its ground state level by emitting a 1809 keV photon. 

The synthesis of ^^Al requires moderate peak temperatures, of the order of Tpeak < 2 x 
10^ K, and a fast decline from maximum temperatures, conditions that are achieved in typ- 
ical nova outbursts |31l43l22j . "^^kX synthesis proceeds through 2^Mg(p, -ff^ k\{l3+ f^Mg{Y>, 
7)2^AF, whereas it is mainly destroyed by (p,7) reactions |22I- A significant nuclear un- 
certainty affects the 2^Al(p,7)2^Si rate |7|22j . which translates into an uncertainty in the 
expected contribution of novae to the Galactic ^^Al content. Calculations based on recent 
prescriptions for the composition of ONe white dwarf cores suggest that the contribution 
of classical nova outbursts to the Galactic ^^Al abundance is small (i.e., less than 15%), 
in agreement with the results derived from the COMPTEL/CGRO map of the 1809 keV 
2^A1 emission in the Galaxy (see [ID]), which points towards young progenitors. 
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2. Observational constraints: from nova shells to presolar grains 

The theoretical nucleosynthetic predictions described above can be compared with the 
abundance patterns inferred from observations of ejected nova shells (see |2H43j . and ref- 
erences therein). The comparison yields in general good agreement between models and 
observations: this includes atomic abundance determinations -H, He, C, O, Ne, Na...-, 
as well as a plausible endpoint for nova nucleosynthesis (around Ca), suggesting that the 
thermal history of the explosion (i.e., Tpeak, exposure times...) is reasonably well repro- 
duced by current models. Unfortunately, observations provide only direct information on 
atomic abundances and hence, do not pose severe constraints on the models. 

Indeed, better constraints can be (partially) obtained from the laboratory analysis of 
presolar grains, which yields isotopic abundance ratios. Classical novae are stellar dust 
factories: infrared measurements in a number of recent novae reveal the presence of C- 
rich dust (Novae Aql 1995, V838 Her 1991, PW Vul 1984...), SiC (Novae Aql 1982, V842 
Cen 1986...), hydrocarbons (Novae V842 Cen 1986, V705 Cas 1993...), or SiOs (Novae 
V1370 Aql 1982, V705 Cas 1993). Remarkable examples, such as novae QV Vul 1987, 
exhibit simultaneous formation of all those types of dust (see jT2] for details on dust- 
forming novae). Recently, several characteristic nova signatures have been identified by 
laboratory isotopic measurements in five silicon carbide and two graphite grains isolated 
from the Murchison and Acfer 094 meteorites These tiny spherules, only a few microns 
in size, are characterized by low ^"^C/^^C and ^^N/-^^N ratios, ^^Si excesses and close-to- or 
slightly lower-than-solar ^^Si/^^Si ratios. In some cases, high ^^A1/^''A1 and low ^''Ne/^^Ne 
ratios have been also determined. This discovery represents the first direct evidence of 
grains that exhibit nova signatures and opens up interesting possibilities for the future. 

Theoretical efforts to predict the expected imprints of nova outbursts on presolar grains 
have been conducted by different authors |42|23j , including preliminary estimates on equi- 
librium condensation sequences in the ejected shells |23j: these studies suggest that clas- 
sical novae may contribute to the known presolar corundum (AI2O3), spinel (MgAl204), 
enstatite (MgSiOs), silicon carbide (SiC) and silicon nitride (Si3N4) grain populations. 

3. Multidimensional effects in nova outbursts 

The assumption of spherical symmetry in nova models (and in general, in stellar ex- 
plosions) excludes an entire sequence of events associated with the way in which a TNR 
initiates (presumably as a point-source ignition) and propagates. The first study of local- 
ized TNRs on white dwarfs was carried out by Shara |HH1 on the basis of semianalytical 
models. He suggested that heat transport was too inefficient for a ffame to spread a 
localized TNR to the rest of the white dwarf surface (i.e., the diffusively propagated 
burning wave may require tens of years to extend throughout the whole stellar surface). 
Therefore, he concluded that localized, volcanic-like TNRs were likely to occur (mainly 
in My^d > I.2M0 white dwarfs). But his analysis, based only on radiative and conduc- 
tive transport, ignored the major (and crucial) role played by convection on the lateral 
thermalization of a TNR. 

The importance of multidimensional effects for TNRs in thin stellar shells and surface 
layers (including classical nova outbursts) was revisited by Fryxell & Woosley |llj . They 
concluded that the most likely situation in nova outbursts involves TNRs propagated by 
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small-scale turbulences. From dimensional analysis and flame theory, the authors derived 
a relation for the velocity of the deflagration front spreading around the stellar surface: 
Vdef ~ (hpVconv / TburnY^'^ , whcrc hp IS the prcssure scale height and Vconv the characteristic 
convective velocity. Typical values for nova outbursts yield Vdej ~ 10^ cm s~^ (that is, 
the flame propagates halfway through the stellar surface in about ~ 1.3 days). 

The first, pioneering studies that addressed this issue in the framework of multidimen- 
sional hydro calculations were performed by Shankar, Arnett & Fryxell [35j, and Shankar 
& Arnett [33]. For that purpose, an accreting 1.25 Mq white dwarf was evolved with a 1-D 
hydro code and mapped into a 2-D domain (a spherical-polar grid of 25x60 km). The ex- 
plosive event was then followed with a 2-D version of the Eulerian code PROMETHEUS. 
A 12 isotope network, ranging from H to ^^F, was included to account for the energetics 
of the explosion. Unfortunately, the subsonic nature of the problem, coupled with the use 
of an explicit code (with a timestep limited by the Courant-Friedrichs-Levy condition), 
posed severe limitations on the study, which was restricted to very extreme (rare) cases, 
characterized by huge T perturbations of about ~ 100 — 600%, in small local regions of 
the envelope's base. The overall computed time was also extremely small (about 1 sec- 
ond). The calculations revealed that instantaneous, local temperature fiuctuations cause 
Ray leigh- Taylor instabilities. The rapid rise and subsequent expansion (in a dynamical 
timescale) cools down the hot material and halts the lateral spread of the burning front, 
suggesting that such local temperature fiuctuations do not play a relevant role in the 
initiation of the TNR (in particular, at early stages). The study, therefore, favors the 
occurrence of the local volcanic-like TNRs argued by 

Later on, Glasner & Livne and Glasner, Livne & Truran ^3], extended these early 
attempts. New 2-D simulations were performed with the code VULCAN, an arbitrarily 
Lagrangian Eulerian (ALE) hydrocode with capability to handle both explicit and implicit 
steps. As in j35l34j . a slice of the star (0.1 tt^"''^), in spherical-polar coordinates with 
refiecting boundary conditions, was adopted. The resolution near the envelope's base 
was around 5x5 km. As before, the evolution of an accreting 1 Mq CO white dwarf was 
initially followed by means of a 1-D hydro code (to overcome the early, computationally 
challenging phases of the TNR), and then mapped into a 2-D domain as soon as the 
temperature at the envelope's base reached ~ 10^ K. As in j3B|34j . the 2-D runs 
included a 12 isotope network. The simulations revealed a good agreement with the gross 
picture described by 1-D models (for instance, the critical role played by the /3+-unstable 
nuclei ^^N, ^'''^^O, and ^''F, in the ejection stage, and consequently, the presence of large 
amounts of ^^C, ^^N and ^"^O in the ejecta). However, some remarkable differences were 
also found: first, the TNR was initiated by a myriad of irregular, localized eruptions at the 
envelope's base caused by convection-driven temperature fiuctuations. Hence, combustion 
proceeds as a chain of many localized fiames (not as a thin front), each surviving only a 
few seconds. Nevertheless, they concluded that turbulent diffusion efficiently dissipates 
any local burning around the core. As a result, they suggest that the fast stage of the 
TNR cannot be localized and therefore, the runaway must spread through the entire 
envelope. Second, contrary to 1-D models, the core-envelope interface is now convectively 
unstable, providing a source for the envelope's metallicity enhancement through a Kelvin- 
Helmholtz instability (a mechanism that bears a clear resemblance to the convective 
overshooting proposed by Woosley |33|)- The efficient dredge- up of CO material from 
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the outermost white dwarf layers accounts for a ~ 30% metal enrichment in the envelope 
(the accreted envelope was assumed to be solar-like, without any arbitrary pre-enrichment 
prescription), in agreement with the inferred metallicites in the nova shells ejected from 
CO novae. And third, larger convective eddies, extending up to 2/3th of the envelope's 
height, with characteristic velocities of Vconv ~ 10^ cm s^^, were found in these 2-D 
simulations. Nevertheless, and despite of these differences, the expansion and progress 
of the TNR towards the outer envelope was almost spherically symmetric (although the 
initial burning process was not). 

Results from other 2-D simulations were published, shortly after, by Kercek, Hillebrandt 
& Truran [25, with the aim to confirm the general features reported in |13|14j . In this 
case, a version of the Eulerian PROMETHEUS code was used. A similar domain (a box 
of about 1800 x 1100 km) was adopted despite a cartesian, plane-parallel geometry to 
allow the use of periodic boundary conditions, was chosen. Two resolution runs, one with 
a coarser 5x5 km grid (as in [T^) and a second with a finer 1x1 km grid, were performed. 
Calculations used the same initial model than in Glasner et al. [14J and yield qualitatively 
similar results but somewhat less violent outbursts (i.e., with longer TNRs and lower Tp^ak 
and fejec), caused by large differences in the convective flow patterns: whereas in [T3j, a 
few, large convective eddies dominated the flow, most of the early TNR was now governed 
by small, very stable eddies (with l^ax ~ 200 km) and, accordingly, more limited dredge- 
up and mixing episodes than in [T3j were found. The authors attributed such differences 
to the different geometry and, more significantly, the boundary conditions adopted in the 
simulations. 

The only 3-D nova simulation to date has also been performed by Kercek, Hillebrandt 
& Truran [2^1 • The run, that adopted a computational domain of 1800x1800x1000 km, 
with a resolution of 8x8x8 km, revealed flow patterns dramatically different from those 
found in the 2-D simulations (much more erratic in the 3-D case): mixing by turbulent 
motions took place on very small scales (not fully resolved with the adopted resolution); 
peak temperatures achieved were slightly lower than in the 2-D case (a consequence of 
the slower and more limited dredge-up of core material). Moreover, the envelope attained 
a maximum velocity that was a factor ~ 100 smaller than the escape velocity and, pre- 
sumably, no mass ejection was expected (except for a possible wind mass- loss phase). In 
view of these results, the authors concluded that CO mixing must take place prior to the 
TNR, in contrast with the main results reported by Glasner et al. |14j . 

4. To explote or not to explote: discussion and outlook 

Despite the thermonuclear runaway model reproduces the gross observational features 
of classical nova outbursts, much remains to be done. The following list (by no means 
exhaustive) outlines different aspects that require further attention: 

- Identify the main mechanism (or mechanisms) responsible for mixing at the core/envelope 
interface. 

- Reanalysis of the expected amount of material ejected per event, as compared with 
the values inferred from observations. 

- Better spectra and analysis technics to provide theoreticians with reliable abundance 
patterns in nova shells. 
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- The observation of ''Li in the ejecta accompanying a nova outburst has been extraor- 
dinarily challenging. Recently, the presence of this elusive isotope has been claimed for 
the first time: an observed feature compatible with the doublet at 6708 A of Li I has 
been reported in the spectra of V382 Vel (Nova Velorum 1999) 9J. However, its has 
been argued that such observed feature in V382 Vel may correspond instead to another 
low-ionization emission centered at around 6705 A likely the doublet associated with N 
I 122] • Confirmation of the presence of this isotope in other nova shells would be highly 
desirable. 

- 7-ray observations of specific gamma-ray signals: that would confirm another of the 
long-standing predictions of the thermonuclear runaway model, either by detecting line 
(478, 511, and 1275 keV) or continuum emission. 

- More theoretical and experimental work required to shed light into the contribution 
of novae to the Galactic ^^Al content (including studies of nova frequencies and galactic 
distributions) 

- New nuclear physics experiments to reduce the uncertainties associated to key reac- 
tions (in particular, ^^Al(p,7) or ^°P(p,7)). 

- Identification of more nova grain candidates. In particular, it would be interesting to 
measure sulfur isotopic ratios (never achieved in presolar SiC grains) and to identify as 
well nova candidates in the oxide grain population. 

- More efforts in numerical modeling, specifically in the multidimensional framework. 
The limited (and somewhat contradictory) efforts to date have to be extended taking 
advantage of state-of-the-art, massive parallel architectures. First, it is crucial to under- 
stand the reason for the differences reported by Glasner et al. and Kercek et al. |21] : 
critical aspects, such as the specific geometry and the boundary conditions adopted, as 
well as numerical inaccuracies in the mapping procedure from 1-D models to a multi-D 
frame can certainly influence the outcome. Moreover, the different flow patterns found 
when comparing 2- and 3-D simulations suggest that probably only 3-D modeling will 
bring final answers to our quest for the real nature of nova outbursts. 
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